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ABSTRACT 

We present time-dependent models of the remnant accretion disks created during 
compact object mergers, focusing on the energy available from accretion at late times 
and the composition of the disk and its outflows. We calculate the dynamics near the 
outer edge of the disk, which contains the majority of the disk's mass and determines 
the accretion rate onto the central black hole. This treatment allows us to follow the 
evolution over much longer timescales (100 s or longer) than current hydrodynamic 
simulations. At late times the disk becomes advective and its properties asymptote 
to self-similar solutions with an accretion rate Md oc (neglecting outflows). 

This late-time accretion can in principle provide sufficient energy to power the late- 
time activity observed by Swift from some short-duration gamma-ray bursts. However, 
because outflows during the advective phase unbind the majority of the remaining 
mass, it is difficult for the remnant disk alone to produce significant accretion power 
well beyond the onset of the advective phase. Unless the viscosity is quite low {a ^ 
10~^), this occurs before the start of observed flaring at ~ 30 s; continued mass 
inflow at late times thus appears required to explain the late-time activity from short- 
duration gamma-ray bursts. We show that the composition of the disk freezes-out when 
the disk is relatively neutron rich (electron fraction Yg — 0.3). Roughly 1O~^M0 of 
this neutron-rich material is ejected by winds at late times. During earlier, neutrino- 
cooled phases of accretion, neutrino irradiation of the disk produces a wind with 
Ye ~ 0.5, which synthesizes at most ^ 1Q~^ Mq of ^^Ni. We highlight what conditions 
are favorable for ^^Ni production and predict, in the best cases, optical and infrared 
transients peaking ~ 0.5 — 2 days after the burst, with fluxes a factor of ~ 10 below 
the current observational limits. 

Key words: accretion disks — black hole physics — gamma rays: bursts — neutrinos 



1 INTRODUCTION 

The most popular model for the creation of short du- 
ration gamma-ray bursts (GRBs) is either binary neu- 
tron star (N S/NS) or black h ole-neutron star (BH/NS) 
coalescence llPaczviiskil 1 19861 . Il99ll : lEichler et aP 1 19891 : 



iNaravan et al.l 19911 '). Support for the merger hypothe- 
sis comes from their durations of 2 s, observations 
of well-localized short GRBs in ga l axies without stroni 
star formation (iBerger et al.l l2005l : iGehrels et "aTl 1200 
iHiorth et al] |2005| ). and the lack of a detectable coin- 
cident supernovae (IHiorth et al.1 20051: Bloom et al.l l2006l : 
ISoderberg et al ] |2006l : lFerrero et al.|[2007l ). as is found in the 



case of lo ng (>, 2 s) GRB s l|Galama et al.lll998l : IHiorth et al.1 
[2003; Sta nek et al.ll2003l ). 

Previous theoretical studies of the merger process have 
focused on one of two stages. The first is the dynamical por- 
tion in which the less mas sive companion is tidally disrupted 
by the more massive BH llLee fc Kluzn iak 1995 , 'l998 l.ll999l: 
Kluzniak fc Lee 19981: Janka et al. _ 1999: Rosswog et al.l 



20041) or NS (iRuffert et al.ril996l : iRuffert fc Jankal Il999l : 



Oechslin fc Jan ka '2006*) . The details o f whether a dyn am- 
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ical instability ( Rasio fc Shapiro! 1 19941 : iLai et al.1 1 1994) ') or 
Roche lobe overflow occurs depe nds on the mass ratio 
and the nuclear equat ion of state (|Bildsten fc Cutierlll992l : 
lUrvii fc Eriguchilll999^ . 

Nevertheless, generally ~ 0.01 — O.IMq of material 
remains in a remnant disk following the dynamical stage. 
The accretion of this material onto the central object gives 
rise to the second, disk portion of the merger. The ener- 
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getics and timescale of the accretion phase are reasonably 
consistent with observations of short GRBs, as was shown 
by models of steady- state, azimuthally symmetric, verti- 



by models oi steady- state, azimutnally symmetric, verti- 
cally a veraged disks JPopham et al. 19991: Naravan et alj 



2OOII: iKohri fc Mincshigc' |2002| : iDiMatteo et alj |2002| : 
Chen fc Beloborodov .2007. '). More recently, these disks have 
been modeled wit h tim e-dependent ca l culation s in ID 
Janiuk et al] 1 2004'). 2D fcee et al.|[20oi . l2005bl ). and 3D 



Setiawan et a l. 2004, 2006). The typical time interval that 



present multi-dimensional calculations can simulate is on 
the order of the burst duration or less ('^ 1 — 2 s for 2D and 
~ 50 ms for 3D). 



Recent observations of short GRBs by Swift, how- 
ever, indicate continued activity from the central engine 
on much longer timescales. X-ray flares with durations 
of ^ 100 s after a delay of ~ 30 s have been seen 
from several bursts (iBarthelmv et al.ll2005l: | Villasenor et al.l 



I2OO5I : ICampana et al.ll2006l : iLa Parola et al.ll2006l) . Stacked 
lightcurves of many bursts indicate cont i nued activity on 
a sim ilar timescale (|Lazzati et al.l I2OOII : iMontanari et al] 
[iooi). In one extreme case, GRB 050724 displayed an X- 
ray flare 12 hours post-burst. This flaring activity has been 
attributed to a number of different sources, in cluding frag- 
mentation of a rapidly rotating core (iKing et ^l. ^2005'1. mag - 
netic regulation of the accretion flow (jProga fc Zhanell200q ) , 
fragmentation of the accretion disk (Perna et al. 2005; al- 
though this explanation may have difficulty reproducing the 
observed timescales, Piro & Pfahl 20 07), di fferential rotation 
in a post- merger millisecond pulsar (iDai et a l. 2006). and an 
infalling tidal tail of material stripped from the disrupted NS 
(Lee & Ramirez-Ruiz 2007; Rosswog 2007). 

In order to determine whether the late-time activity 
from short GRBs is consistent with a compact merger origin, 
the disk evolution should be followed for timescales much 
longer than the initial viscous time. With this aim, we per- 
form time-dependent calculations modeling the disk as an 
annulus that contains the majority of the mass. This simpli- 
fication allows us to study the disk evolution for arbitrarily 
long timescales, and to readily determine important prop- 
erties such as the disk's composition and when it becomes 
advective. We are also able to survey much of the param- 
eter space of initial disk mass and angular momentum. In 
32] we discuss the initial conditions for disks formed from 
compact object mergers. This is followed by SjSl in which we 
summarize the main assumptions of our ring model. In SJJ] 
we present the results of our calculations and summarize the 
main properties of the models. We then calculate outflows 
from our disk solutions in SJS] We investigate the compo- 
sition of the outflows and argue that they generally con- 
sist of neutron-rich isotopes, but can produce ^®Ni in some 
circumstances. The presence or lack of an optical transient 
from short GRBs therefore provides an important constraint 
on progenitor models. We conclude in JHlwith a discussion 
of our results. In Appendix A we summarize the Green's 
function solution to the viscous spreading of a ring, which 
is important for connecting our ring model to the true ex- 
tended disk geometry. In Appendix B we present analytic 
self-similar solutions that reproduce many of the features of 
our numerical solutions. 



2 INITIAL CONDITIONS 

The dynamical phase of NS/NS or BH/NS mergers has 
been studied extensively using a number of different nu- 
merical techniques and methods for including general rel- 
ativity (GR). Here we summarize some of the most rele- 
vant f eatures for our study (for a more detailed review, see 
iLee fc Ramire z-Ruiz 2007i). 

When the lighter companion NS is first tidally dis- 
rupted, a debris disk is formed within only a few dynamical 
timescales. The initial disk mass, Md,o, is generally larger 
for more asymmetric mass ratios (i.e., small q, where q is 
the ratio of t he lighter to the heav i er bin ary component). 
For example, IShibata fc Taniguchil l|2006l ) find that for a 
NS/NS merger with q = 0.7 that Md,o = 0.03 M©, but for 
q — 0.9 the disk is much less massive with Md,o = 10~^ Mq. 
Another trend is that including strong gravity gives less 
massive remnant disks. The BH spin i s also important , 
with larger spin favoring disk formation (|Rasio et alJbOOSh 
and the production of a tidal tail. These have masses of 
^ 0.01 — 0.05 Mf7) and may provide prolonged mass inflow 
jLee fc Ramirez-Ruizir2007l ). but for simplicity this will be 
ignored here. Taken together, these simulations generally 
flnd Mdfi — 0.01 — O.3M0, with the disk containing a sub- 
stantial fraction of the angular momentum of the disrupted 
companion. 

In the standard picture of NS-NS mergers, the result- 
ing hypermassive NS collapses to a BH shortly following the 
merger. However, simulations show that when (and if) col- 
lapse actually occurs depends on the mass of the central 
NS and its ability to transport angular momentum to the 
surrounding disk (Shibata et al. 2005; Shibata fc Taniguchi 
2006; Shibata et al. 2006). In fact, if the NS remains sup- 
ported by differential rotaton for several seconds (Baum- 
garte et al. 2000; Morrison et al. 2004; Duez et al. 2004, 
2006) or loses sufficient mass via a centrifugally-driven out- 
flow (e.g., Thompson et al. 2004; Dessart et al. 2008a), the 
NS itself may power the GRB (e.g.. Price fc Rosswog 2006). 
In this paper we assume that the central object promptly 
collapses to a BH; our model, however, would be reasonably 
applicable for the case of a central NS as well, the primary 
difference being that the significant neutrino fiux from the 
newly-formed NS and from the boundary layer between the 
disk and the NS could modify the composition and thermal 
properties of the disk. 

We present some characteristic numbers to motivate our 
choice of initial conditions. Consider a binary with masses 
Al and m (M > m), where the latter is the NS (with radius 
R) that is tidally disrupted. The disruption radius, at, is 
estimated to be (Kopal 1959, adding Fishbone's 1973 10% 
strong gravity correction) 



at ~ 2.47? 

V m / 

The characteristic orbital period at this radius is 

,3 \ 1/2 



Pt ~ 23.4 



2 X 10~^mi y^i?g''^ s, 



(1) 



(2) 



where mi. 4 = m/1.4 Mq and i?6 
orbital angular momentum of 

Jt = {G{M + m)atf'^m 



i?/10^ cm, with an 
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ergs s, 



(3) 



3.1 Dynamical Equations 



where q — m/M. The disrupted NS also contains spin angu- 
lar momentum. This is neg ligible since the NS is not strongly 
affected by tidal coupling (|Bildsten fc Cutlerlll992l '). Even a 
rapidly rotating NS (~ 5 ms) has an associated angular mo- 
mentum of merely ~ 10** ergs s. 

Once disrupted, a considerable fraction of the NS is ei- 
ther lost from the system or immediately swallowed by the 
BH. The remaining material forms a thick torus surround- 
ing the central BH. Its associated viscous timescale can be 
estimated by assuming that the majority of the torus' mass 
lies at a single radius, rd,o- Taking the angular momentum 
of the disk to be Jd — {GMrdfiY^'^ Md,o, we estimate 

r,,o ~ 3 X loV/g-^MoTi" ' cm, (4) 

where M3 = M/3 Mq, Mq.i = Mdfi/O.l Mq, and J49 = 
ergs s. For a disk with half-thickness H, the viscous 
timescale is 



^visc.O 



6 X 



\hj \gm ) 

in-2 -l,,r-l/2 3/2 ( H \ 

1° los^j 



(5) 



where a = O.lap.i is the standard dimensionless viscosity 
IShakura fc Sunvaevlll973l ). = r^.o/lO^ cm, and we have 
scaled to an initial ratio of H/vd = 0.5, consistent with 
our numerical solutions. The initial viscous time tvisc.o is 
roughly the time at which the central BH begins accreting 
in earnest. The strong dependence of fvisc.o on disk mass 
and radius demonstrates that the initial evolution of the 
disk is sensitive to the outcome of the dynamical phase of 
the merger. But as we will show, the late time evolution is 
much less sensitive to initial conditions and is well described 
by self-similar solutions. 



3 PHYSICS OF THE EXPANDING RING 
MODEL 

Given these initial conditions, one would like to know how 
the disk then evolves. Modeling the entire disk requires re- 
solving timescales over ~ 4 — 6 orders of magnitude. This 
makes it expensive to carry out simulations for long periods 
of time. We consider instead a simplified model that cap- 
tures most of the features of interest. At any given time, 
t, the disk can be broken into three regions depending on 
the local viscous time, tvisc, which increases with radius, 
roughly as tvisc ~ r^^^. At small radii, fvisc < t, and the disk 
comes into steady-state. This is the region most often mod- 
eled i n previous studies |Popham et al. 19991: Naravan et al.l 



200ll: iKohri fc Mineshige 
Chen fc Beloborodovl 12007 



2OO2I : iDiMatteo et al.l booa 

The radii where fvisc ~ t con- 
tain the majority of the disk's mass and angular momentum. 
Therefore, this region determines the viscous evolution of 
the rest of the disk, including the mass accretion rate that 
is fed to the interior steady-state region. Motivated by this 
fact, we focus on this radius and model the disk as a ring. 
Exterior to this point is a third region where fvisc > t, but 
this contains a small amount of mass and is negligible for 
the viscous evolution. 



Our ring model treats the disk as a single annulus that is 
evolved forward in time. In this picture, the properties of the 
ring, such as its surface density S and temperature T, are 
representative of the location where Er^ peaks. The main 
drawback of this method is that the material in the disk is 
in fact distributed spatially in radius. Thus, although the 
mass of the disk in the vicinity of is ~ vrEr^, the total 
mass of the disk (integrated over all radii) is Md = AyrEr^, 
where A is a factor of order unity that accounts for the dis- 
tinction between the total mass of the disk and the mass 
of the material near r^. Similarly, we write the total angu- 
lar momentum of the disk as Jd — B{GMrd)^^^'!^rd^. At 
early times the constants A and B depend on the initial 
conditions of how matter is spatially distributed; however, 
at times much greater than the initial viscous time (given by 
eq. [S]), material initially concentrated at a given radius be- 
comes spread out in a manner determined by the viscosity. 
As described in detail in Appendix A, we choose the con- 
stants A and B by setting the solution of our simplified ring 
model at late times equal to the Green's function solution 
for a spreading ring with a viscosity 1/ cx r^^'^ (as is appro- 
priate for the radiatively inefficient disk at late-times). This 
fixes A = 3.62 and B = 3.240 Conveniently A/B ~ 1, so 
that it is a good approximation to take Jd — {GMrd)^^^ Md- 
The time evolution of the disk is determined by the 
conservation equations. Conservation of mass is 



dt 



[ATiEr^d] = 



'Md 



(6) 



where Md is in general the total mass loss rate, which could 
include both accretion and a wind (for now we ignore the 
effects of a wind) . Conservation of angular momentum is 



£ [B{GMrdy^\^rl] = - j. 



dt 



(7) 



where J is the angular momentum loss rate. Equations ([6} 
and (O provide two coupled equations that can be solved 
for the dependent variables and E. 

The accretion rate must depend on the characteristic 
mass and viscous timescale of the ring, so we use 



Md = fMd/Ui,: 



(8) 



where tvisc ~ r'^/v and v is the viscosity. The factor / is 
set like A and B to match the exact solution of a spreading 
ring with u cx r^''^ (Appendix A), which gives / — I.gQ 
Requiring a no-torque boundary condition at a radius r«, 
we take 



f=l.Q/[l~{r./rdf'% 



(9) 



In contrast, a steady-state disk obeys Md ~ SnvT, (ignoring 
the no-torque condition), which instead gives / = 3/ A ~ 
0.83. 

For the viscosity, we use an a-prescription, 



u = aCsH, 



(10) 



^ In fact, when the total angular momentum is conserved, the 
viscous evolution is independent o{ A/B as long as A/B is nearly 
constant with time. 

^ Although we set tvisc = /^, any prefactors that could go into 
this prescription would just be absorbed into a re-definition of /. 
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where Cs = {P/pY^"^ is the isothermal sound speed. The 
equation of state includes contributions from radiation pres- 
sure, gas pressure, r elativistic degene racy pressure, and neu- 
trino pressure as in iDiMatteo etal] (fiooa l. 



3.2 Energetics 

For the energy equation, we take 

5visc ~ ?,y + Qadv , 



(11) 



where (/vise is the viscous heating, q' J is the neutrino cool- 
ing (using the prescriptions given by IDiMatteo et al]l2002l . 
which includes neutrino optical-depth effects), Qadv is the 
advective heat flux, and all q values correspond to half the 
disk thickness. 

For a disk rotating at the Keplerian frequency Q. — 
{GMIrlY'\ 



9 9 GMMd 

8 8/4 7rr3 



1/2 



(12) 



where the prefactor 9/(8/yl) ~ 0.2 is different from the 
stead y-state value of 3/8. T he advective term, gadv, is set 
as in IDiMatteo et al] (|2002f ). with the only difference being 
that the radial velocity is the expansion rate of the ring's 
radius 



Vr 



dt 



2M 



(13) 



where we have taken J = 0. 

Fusion to Qf-particles produces heating in addition to 
gvisc, with 

gnuci = 6.8x lO'Vio^^^, (14) 

where all quantities are expressed in cgs units, pio = 
p/10^° g cm~^ and Xa is the mass fraction of a-particles. 
Note that in our case Qnuci > since a-particles are synthe- 
sized as the disk expands (in contrast to studies that follow 
cooling from photodisintegration as material moves inward) . 
In our calculations we do not include gnuci in solving equa- 
tion (|11|) because we were not able to find reasonable solu- 
tions when doing so (for reasons explained in tj4.ip . 



3.3 Composition 

An advantage of the ring model is that other properties of 
the disk, such as its composition, can be cast into differential 
equations and integrated along with equations (|6} and ([7]). 
Since the neutron content of the disk is particularly impor- 
tant for determining the properties of the disk's outfiows, 
we evolve the electron fraction using 



dt 



(15) 



where = X^j (X„ 4- Xp), X^ and X^ are the proton and 
neutron mass fraction, respectively, and r^-^ and are 
the electron and positron capture rates, respectively (Be- 
loborodov 2003a). We have neglected the effect of neutrino 
absorptions on the evolution of Vg in equation (|15p . Al- 
though absorptions are important at early times when the 
disk is optically thick, we are primarily concerned with the 



late-time value of Fe, which does not depend sensitively on 
the neutrino irradiation (see >j4.2|l . 

As the disk evolves, the protons and neutrons eventu- 
ally burn to form a-particles. At these times the disk is 
sufficiently cold that the positron and electron capture rates 
are negligible (i.e, l/r^-p 3> fvisc) and Fe has frozen-out. 
This fixes the difference between the free neutron and pro- 
ton mass fractions 



Xn — X^ — 1 — 2ye. 



(16) 



Since the rates for reactions that synthesize and destroy 
a-particles are all fast in comparison to the viscous time, 
we determine the composition using nuclear statistical equi- 
librium (NSE) between protons, neu tron, and a-particles. 
This is expressed by the Saha relation jShapiro fc Teukolskvl 
119831 ) 

Xlxl = 1.57 X 10* X^pl^Tli'' exp (-^) • (17) 

NSE is a good assumption because the disk temperature 
is generally > 0.5 MeV (see Fig. [5}, except at very late 
times or for very low disk masses (e.g., the Md,o = 0.03Mq 
case, for which we do not calculate the nuclear composition 
anyways). By combining equations (|16|l and (|17p with mass 
conservation, Xp -fX^ = 1, we solve for all of the mass 
fractions at a given p, T, and Ye- 



4 TIME-EVOLVING SOLUTIONS 

We next present the results of integrating equations (O, (O, 
and (|15p forward in time. For simplicity, we typically assume 
that J = 0. A convenient property of our formalism is the 
ease with which these complications can be included (for ex- 
ample, we consider the effects of winds at the end of §4.1). 
The disk properties are determined by the initial conditions 
Md,Q, Jd, and Ye,o, and by the viscosity a. For the majority 
of our study we set the initial Y^fi = 0. 1, which is charac- 
teristi c of the inner neutron star cr ust (|Haensel fc Zduiii^ 
Il990al lbl: |p"ethick fc RavenhaU Ill995l l. An additional impor- 
tant parameter is r*, which is set by the spin of the cen- 
tral BH. In most of our calculations we take r* ~ 2.3r5 ~ 
1.02 X 10*' cm, corresponding to the innermost stable circu- 
lar orbit of a 3 Mq BH with spin a ~ 0.9; when calculating 
the properties of disk outfiows in ij5l however, we also con- 
sider the case of a nonrotating (a = 0) BH. We consider the 
general evolution of the disk in ^4.11 and then focus on the 
composition in M.2\ 



4.1 Disk Evolution and Energetics 

At any given time, a ring model is in one of three phases: 
(1) early-time, optically thick to neutrinos and advectively 
dominated, (2) mid-time, optically thin to neutrinos and 
geometrically thin, and (3) late-time, radiatively-inefficient 
accretion flow (RIAF)I!| This is analogous to the different 
regions of steady-state, hyper-accreting accretion disks (see. 



An optically thick, geometrically thin stage occurs between 
stages (1) and (2); however, this phase is brief and is not dy- 
namically very different from phase (2), so we do not consider it 
separately in our discussion. 
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Figure 1. Example disk models siiowing the evolution of the disk 
radius, r^, disk mass, Md, and accretion rate, M^, as a function 
of time. We compare M^j g = 0.03 (solid lines), 0.1 {dotted lines) 
and 0.3 Mq (dashed lines) solutions; all use J49 = 2 and a = 0.1. 
The inner radius is r* ~ 2.3rg ~ 1.02 X 10*5 cm (corresponding to 
a 3 Mq BH with a spin of a ~ 0.9). 




0.2 



-2 
logio t (s) 



Figure 2. Comparison of the midplane temperatures and scale- 
height for the three models from Fig. 1. In the lowest mass model, 
the ring is always advectively-dominated, thus H/r^ is constant. 



Figure 3. The cooling rates and neutrino luminosity for the 
Md.o = O.3A^0 model from Fig.[T] For the cooling rates we com- 
pare the neutrino {solid line) and advective {dashed line) rates, 
normalized to the viscous heating. The implied heating from the 
creation of a-particles is plotted as a dotted line, but is not ac- 
counted for in the disk evolution. The neutrino luminosities are 
from the entire disk {solid line) and the ring {dashed line). The 
former luminosity is estimated by integrating over a steady-state 
disk model at each time given Mrf(t). 



e.g., IChen fc Beloborodovl |2007| ). but now the transitions 
occur with time instead of radius. The phases that a cer- 
tain ring model samples during the course of its viscous ex- 
pansion depends on tvisc.o- A more compact disk (a shorter 
tvisc.o) will exhibit all three phases, while larger disks may 
only exhibit phases (2) and (3), or even just (3). 

We present a number of figures that are helpful in un- 
derstanding these three phases and how they are affected by 
changing M^fi- Figure [T] shows the radius r^, mass Md, and 
accretion rate Add as a function of time, for Mdfi ~ 0.3, 0.1, 
and 0.03 A/0. Figure [2] compares the midplane temperature 
and scaleheight for these same models. Figures [3] and [3] show 
key results describing the energetics of the Md,o = 0.3 and 

0. 1 Mq solutions, respectively, while Figure [5] shows the dif- 
ferent contributions to the total pressure in the disk as a 
function of time. Note that we fix the total angular momen- 
tum in these calculations (J4g — 2) and thus a larger A'Id,o 
corresponds to a smaller rd,o and a shorter fvisc.o- 

The first transition the disks make is from an opti- 
cally thick, advective disk to a thin, neutrino-cooled disk; 

1. e., from phase (1) to (2). This is only exhibited by the 
Mdfi = O.3M0 model and is seen most clearly at early 
times in Figure [5] when H/vd — 0.5 and in Figure [3] when 
gadv ^ ■ Figure [5] shows that this phase is ion pressure 
(ideal gas) dominated. A simple estimate determines what 
initial disk mass is required for phase (1) to occur, i.e., for 
the initial disk to be optically thick and advective. The disk 
is advective for radii inside of which the neutrino diffusion 
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logio t (s) 

Figure 4. The same as Fig. [3] but for M^fi = 0-IMq. 

time out of the disk exceeds the inflow time. Setting this 
radius equal to the initial radius of the disk (eq. [4]), we find 
that there is a critical disk mass below which the disk never 
experiences phase (1), 

M.,„. ~0.2a-/-M3-/- (-|-)"''Me,(18) 

where we have dropped scalings with / and A since they 
appear raised to the 1/10 power. This estimate is consistent 
with the fact that our Md,i = O.IMq model is not advective 
at early times, as seen in Figures [2] and |4] In this case only 
phases (2) and (3) are seen, i.e., the disk is initially thin and 
neutrino cooled and later transitions to being advective. 

Once the models reach the late-time, RIAF phase, or 
phase (3), they asymptote to self-similar solutions, indepen- 
dent of the initial disk mass. In this phase, the disk has 
^Zadv > <lZ a-iid is radiation pressure dominated. We derive 
analytic self-similar solutions in Appendix B2 for this limit 
and show that cx i^^^, Md cx t~^l^ and Md oc t""/^. The 
RIAF solution occurs external to an "ignition radius," which 
we estimate as the location where the pair capture cooling 
rate balances ~ 1/2 of the viscous heating for a thick disk, 

^ign — 

/ , X -14/5 / ,v \ 6/5 

3xl0ao,M3 / (^^j (^^^YMei^j ^^'^ 

where we have scaled H/r^ to ~ 0.4, a value appropriate for 
the transition between the thin and thick disk regimes. We 
combine this with the analytic results for rd{t) and Md{t) 
in the RIAF limit (eqs. |B7j and [BBI 'FI to estimate the time 
when the disk transitions to being thick, which yields 

* We use these solutions rather than the thin-disk ones because 
the numerical results follow these more closely (Fig.lBTIl. 
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Figure 5. Pressure contributions for o = 0.3Mq (top panel), 
0.1A/q {middle panel) and O.OSMq {bottom panel). The pressures 
are all normalized to the total pressure and include the ion pres- 
sure {solid lines), radiation pressure {dotted lines), degenerate 
electron pressure {dashed line), and neutrino pressure {dot-dashed 
line) . 

, ni -23/17 ,,-13/17 / J49 ,„„^ 

tthick ~ O.IQq 1 M3 I ~2~J ^' ^ ' 

Equation (|20|) is only applicable if the disk is thin at early 
times. For sufficiently small initial disk masses, less than 

Md,tMck ~ O.lal^fM-'^'' (:^) Mo, (21) 

this is no longer true, and the disk is always a RIAF at its 
outer radius. 

Figures [3] and |4] show that at approximately the same 
time as the disk transitions from being thin to thick, protons 
and neutrons are fused to He. Although the nuclear heating 
rate gnuc is shown in Figures [3] and |4l this heating was not 
included in our time-dependent calculations so that we could 
obtain solutions at late times. The nuclear heating rate is 
sufficiently large, i.e, Qnuci > Qvisc, that the disk is not able 
to accommodate this added energy (it is already thick with 
H ~ r due to viscous heating alone). This probably implies 
that the bur ning contributes to driving a powerful wind (as 
described bv lLee fc Ramirez-Ruij|2007h . 

However, such a wind already begins at this time by 
virtue of the disk being advective (as discussed in i]5.2[) . In 
Appendix B3, we present analytic self-similar solutions for 
advective disks with mass loss and show that this significant 
mass loss causes Md and Md to decline much more rapidly 
with time than is shown in Figure [T] This is shown explic- 
itly in Figure (6] where we present disk models calculated 
using the mass and angular momentum loss prescriptions 
described in Appendix B3; such losses are assumed to oc- 
cur only when the disk is thick, between ^ max(r, ,rign) 
and ~ Vd- Figure |S] compares time-dependent solutions with 
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Figure 6. The radius r^j, disk mass M^, and mass accretion rate 
reaching the central BH Min for different parameterizations of 
mass loss during the advective phase. We initialize a disk with 
^d,o = O.IMq (and all other parameters fixed as in Fig. [TJ and 
compare solutions with no wind [solid line), p = 0.5 {dotted line; 
see eq. [B8]), and p = 1 {dashed line). 

no wind {solid line), a wind with p = 0.5 {dotted line; see 
eq. [B8]), and a wind with p = 1 {dashed line)^ The loss 
of angular momentum does not appreciably slow the radial 
expansion of the disk, but it does substantially accelerate 
the decline in the disk mass and accretion rate (see also eqs. 
[BT2] and [BlS] '). If the models with winds are accurate, sig- 
nificant accretion is only likely to last for a few viscous times 
once the disk enters the late-time advective phase. Contin- 
ued central engine activity at much later times could result 
from late-time infall of ti dally stripped NS material (e.g., 
iLee fc Ramirez-Ruizll2007l ). 

As an additional comparison, we present the effect of 
varying Jd in Figure [T] The main trend is that a higher 
Jd has a larger initial radius for a given Md, and therefore 
a longer viscous time and smaller accretion rate. The late 
time behavior is more sensitive to Jd than the initial Md, as 
predicted by the self-similar solutions, but it still does not 
affect the late time disk radius (see eq. |B7) ). We do not plot 
our results for different a since they are generally consistent 
with the analytic scalings above and in Appendix IB] 

4.2 Composition 

The composition of the disk is important for determining 
the observational effects of any outflows. To this end, we 
plot the composition of our Md,o = 0.3 Mq, J49 — 2 disk 
as a function of time in the upper panel of Figure |H] In 
the bottom panel we plot the relevant timescales for setting 

^ See Appendix IB3I for the definition of p. 
© ???? RAS, MNRAS 000.[TlfT9l 
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Figure 7. Similar to Fig. [T] but now taking the angular momen- 
tum to be J49 = 2 {solid lines), 4 {dotted lines), and 6 {dashed 
lines). All solutions take Md = O.SMq with all other variables 
the same as in Fig. [T] 

the composition, namely the viscous timescale, tvisc {solid 
Itne), the neutronization timescale t„ = l/r^-p {dotted line), 
and the timescale for a-particle photodisintegration, tphoto 
{dashed line). At early times t„ <C tvisc, so that an equilib- 
rium value of Ye ~ 0.23 is reached almost immediately. As 
the disk leaves the optically thick phase and becomes thin- 
ner, degeneracy pressure plays a larger role. This enhances 
neutron production, with a minimum Ye — 0.05. As the neu- 
trino cooling subsides and the disk becomes thick again. Ye 
increases. Before Ye can reach ~ 0.5, it freezes-out at a value 
of Ye ~ 0.3 once t„ > t^iac- 

Besides the neutron abundance. Figure |S] also highlights 
the production of a-particles. Initially, the reactions needed 
to convert neutrons and protons to helium as well as pho- 
todisintegration of helium all happen on timescales much 
shorter than the disk evolution timescale (as an example, 
we plot the helium photodisintegration timescale in the bot- 
tom panel of Fig. so that we can estimate the a-particle 
mass fraction using chemical balance (eq. [IT])- Once the a- 
particle photodisintegration timescale becomes sufficiently 
long (tvisc < tphoto), chemical equilibrium no longer applies 
and Xp = 0, X„ = l- 2Ye ~ 0.4, and = 2Ye ~ 0.6. 

Figure[9]shows how the late-time, frozen-out value of Ye 
in the disk depends on the initial disk mass Md,o and radius 
rd,o, for two different initial electron fractions, Ye,o ~ 0.1 
and Ye,o = 0.5. The former is relevant for the disks cre- 
ated from NS-NS or BH-NS mergers (the focus of this pa- 
per), while a larger Ye,o — 0.5 is appropriate for disks cre- 
ated during the accretion-induced collapse of a white-dwarf 
to a neutron star (e.g., Woosley & Baron 1992; Dessart 
et al. 2006) . Figure [5] shows that for sufficiently compact 
disks, the disk reaches a modestly neutron-rich composition, 
with Ye ~ 0.3 — 0.4, independent of the initial composition. 
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Figure 8. The composition and important reaction timescales 
as a function of time, for the Mdfi = 0.3Mq model from Fig.[T] 
In the top panel we plot the electron fraction, Ye and the mass 
fraction of protons, neutrons, and a-particles {see inset key). In 
the bottom panel we show the viscous time, tvisc {thick, solid 
line), the neutronization time, tn = l/r^-^ {dotted line), and the 
o-particle photodisintegration time, tphoto {dashed line). 

This is because, as highlighted in Figure O the timescale to 
come into /3-equihbrium is shorter than the viscous time. For 
disks with a small initial mass and/or a large initial radius 
(the lower right-hand corner of each panel), t„ > tvisc and 
the disk retains its initial composition (set by the tidally- 
disrupted progenitor and the subsequent dynamical stage of 
the merger). Finally, neutrino irradiation of the outer disk 
by the inner disk can increases the freeze-out electron frac- 
tion, but we estimate this changes the freeze-out value of Ye 
by at most ~ 20%|f| 



5 DISK WINDS 

Having described the evolution of the accretion disk as a 
function of time, we now discuss the properties of outflows 
from these hyper-accreting disks. Winds driven from deep 
within the BH potential well could produce relativistic jets 
and power late-time central engine activity. Outflows driven 
from larger radii dominate the system's mass loss and may 
power supernova-like optical transients through the decay of 
radioactive isotopes that are synthesized in the wind (Li & 
Paczyhski 1998; Kulkarni 2005). In both cases, the mass loss 



° Our calculations employ the pair-capture cooling prescription 
of DiMattco ct al. (2002), which assume Ye = 0.5 and ultra- 
relativistic electrons; we find, however, that including the ef- 
fects of degeneracy and arbitrary electron energies on the cooling 
changes the asymptotic electron fraction by at most a few percent. 




1 10 
^d.o (10^ cm) 

Figure 9. Contours of late-time electron fraction in the expand- 
ing disk as a function of initial disk mass Af^j g a^nd radius r^j q, for 
two different initial compositions. Relatively compact disks come 
into /3-equilibrium and reach an electron fraction independent of 
the initial Ye, while low mass, more extended disks retain their 
initial composition. Figure [S] shows the evolution of Ye with time 
for one particular disk solution. 

rate and nuclear composition are critical for determining the 
observable signature. 

The type and character of the outflow depends on the 
disk's thermodynamic state and changes as it passes through 
the different stages of evolution described in the previous 
section. In i]5.1l we discuss early times when winds are due 
to neutrino irradiation of the thin, efficiently neutrino-cooled 
portions of the disk. We then consider thermally driven 
winds during thick, radiatively-inefficient accretion in ^5.2\ 
This dominates the mass loss at late times and blows away 
most of the remaining disk. In i]5.3l we summarize the nuclear 
composition of the outflows during each phase. We predict 
an ejected '"^Ni mass of at most ~ 10"^Mq {^EM- Its de- 
cay may power transient emission detectable following some 
short GRBs. 

5.1 Neutrino-Heated Thin-Disk Winds 

A wind with a mass loss rate Mm driven from a thin 
disk at radius r must absorb a net power greater than 
Eb = GMMw/2r to become unbound from the central BH. 
In principle, Eb may be supplied by dissipation of the tur- 
bulence that produces the accretion shear stresses. "Vis- 
cous" heating of this kind only efficiently drives an out- 
flow if a substantial fraction of the accretion power is dis- 
sipated in the disk's upper atmosphere, where the cooling 
timescale is long compared to the wind's outward advec- 
tion timescale. However, local radiation MHD simulations 
to date suggest that very little energy dissipation occurs 
in the corona (e.g., Krolik et al. 2007). Instead, heating in 
the atmosphere above a thin, neutrino-cooled disk is likely 
dominated by neutrino irradiation. We therefore focus on 
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the neutrino-driven mass loss rate, which sets a minimum 
Mw, and which can be reliably estimated. Neutrino-driven 
outflows from hyper-accreting disks have also been studied 
by Daigne & Mochkovitch (2002), Levinson (2006), Metzger 
et al. (2008a), and Barzilay & Levinson (2008); Dessart et 
al. (2008b) calculate the neutrino-driven mass loss from the 
central NS following a NS-NS merger under the assumption 
that collapse to a BH is not prompt. 

The neutrino-driven mass loss rate is calculated by 
equating Eh to the total neutrino heating rate in the disk's 
atmosphere. For the radii and entropies that characterize the 
winds, heating via electron neutrino absorption on baryons 
{p+Ve — > n + e^ and n + Ve — >p-|-e~) dominates other forms 
of neutrino heating (e.g., v — u annihilation and v — e~ scat- 
tering; see Qian & Woosley 1996; hereafter QW96). Since 
the neutrino absorption cross section, cti^jv ~ 5 x 10~*'*(eJ) 
MeV~^ cm^, increases with neutrino energy, neutrinos radi- 
ated from near the inner radius r* dominate. Assuming that 
the i^e and Ve. luminosities and spectra are approximately 
equal and can be approximated as originating from a point 
source at small radii, the neutrino heating rate through a 
surface density E at radius r is 



LvOuN^ r,win39r / 2 vv^ -2 -1 

^ TZZ: — I? — 2x10 L52(eio)Ei8r6 ergs s cm 



47rmjvr2 



.(22) 



where r = lOVg cm, = 10''''L52 ergs s-\ (ej) = 100{efo) 
MeV^, and E = EiglO^^ g cm~^. This expression assumes 
that the absorbing layer is optically thin, i.e., that = 
Ea^N/ruN — 3Ei8(eio) < 1- 

First, consider neutrino heating in comparison to vis- 
cous heating in the midplane. This ratio is largest when the 
disk is marginally optically thick (r^ ~ 1), peaking at a 
value of 



Qvisc 



-(e)(t^)(s 



A 



2/5 ,2/5^ 



-6/5 



(23) 



where e = L^/Md(? is the disk's radiative efflciency. Thus, 
although we neglected neutrino heating in S|4l it may become 
somewhat important when ^ 1 and should be included 
in a more detailed calculation. 

We now consider a wind that emerges from the disk 
in the z-direction, parallel to the rotation axis. Away from 
the disk midplane, neutrino heating dominates over viscous 
heating, balancing cooling (g^ = g^J ) at a slightly lower tem- 
perature, ~ S.SLsf (e?o>^''^^-^^''^MeV. Moving further 
out in the hydrostatic atmosphere, the temperature slowly 
decreases below T^. Due to the strong temperature depen- 
dence of the pair capture cooling rate (q~ oc T®), a "gain 
region" of net neutrino heating (i.e., > q^) develops 
above a height Zgain. This net heating drives an outflow. 

The thermal power deposited in the upper disk atmo- 
sphere Eu is the specific heating rate g^/E (eq. [22]) mul- 
tiplied by the mass of the atmosphere in the gain region 

Mgain ~ 27rH(2gain)r^p(%ain), whcrC /^(Zgain) IS the SCalc 

height near the base of the gain region. Although the mid- 
plane of a neutrino-cooled disk is generally dominated by 
nonrelativistic gas pressure (see Fig. [SJ, the gain region has 
a sufficiently low density that it is instead dominated by 
radiation pressure P^ad ~ (ll/12)aT'*. Its scale height is 



^/(zgain) ^ (^'rad/p5z)|zg^i„, where gz is the gravitational ac- 
celeration in the z-direction. Since -ff(2gain) is less than the 
midplane scale height H, Zgain — H and ~ GMH/r^. The 
atmosphere in the gain region is roughly isothermal so we 
set r( Zgain) ~ T^- By combining these estimates and equat- 
ing with Et we find that the neutrino-driven mass loss 
rate from a thin disk is 

Af.|s"»s« ^ 10-'Lli'{elof^\l/'M;^\H/r)-'MQS-\i2A) 

analogous to that derived by QW96 for proto-neutron star 
winds. The assumption that the atmosphere is radiation 
dominated is only valid if the asymptotic entropy in rela- 
tivistic particles exceeds that in nonrelativistic nucleons 
Sn ^ 6 + ln{T^^y/pio)kB baryon-\ where T = TmcV MeV. 
By dividing the energy gained by a nucleon in the wind 
~ GMtun /'2r by the gain region temperature r(zgain), we 
estimate 



5"" ~ QOL7j^''{eloy'-'''r-^^''Ms fca baryon"' 



(25) 



as the asymptotic wind entropy. 

Although equation (|24|) does not strictly hold when 
5"" ^ Sn, QW96 show that AIi, scales the same way with 
L^, (ej), M, and r, but with a larger normalization of 



10" 



5 r 5/3 ,2 X 5/3 5/3, r- 
-'^52 \<:iO/ 



\H/r)-'MQS-\{26) 



The mass loss rate is higher for low entropy winds because 
neutrino heating peaks further off the disk surface, which 
reduces the binding energy and gravitational acceleration of 
matter in the gain region. Using the numerical disk wind 
calculations described in Metzger et al. (2008b; hereafter 
M08b) we have verified that equation (|26[) holds to within a 
factor ~ 2 when S° ~ Sn- 

In deriving equations (|24|) and (|26|l . we have implicitly 
assumed that the timescale for neutrinos to heat matter in 
the gain region thcat = {Uth^/ pqt)\zg.^i^, where Uth 3Prad 
is the thermal energy density, is short compared to ivisc, the 
timescale over which the disk properties appreciably change. 
Equating (eq. [25) ) to the entropy in relativistic particles 
(X T^/p, we find that 

p(zgain) ^ lOV-^/^L^f (e?o)'/'M3-ig cm'^ (27) 
Then, using equations (|22p and (|27p . we have thal[3 



tho 



3-Prad 



P(9.+ /S) 



(28) 



For most of the disk solutions considered in this paper, we 
find that thcat < ivisc during the thin disk phase; thus, equa- 
tions (|24p and (|26|) are reasonably applicable near ra- 

Figure [10] compares the accretion rate Md {solid line) 
with the neutrino-driven mass loss rate M^. In order to de- 
termi ne L^, and (e^), we cal culated steady-state disk models 
(e.g., iDiMatteo etal\ \200± with the accretion rate set at 
each time according to our ring model with J49 = 2 and 
Md = 0.3Mq. We plot the neutrino-driven mass loss rate 
Ml, (eqs. [24] and [26]) at small {dotted line) and large {short- 
dashed line) radii. This shows that the mass loss is dom- 
inated by large radii where the majority of the mass lies. 



^ Equation I I28I I is also approximately equal to the outward ad- 
vection timescale of the wind in the heating region. 
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Figure 10. The accretion rate Af^ {solid line) and neutrino- 
driven mass loss rates Mi, for our J49 = 2 and M^j q = 0.3Mq 
model, focusing on the phase of thin, efficiently neutrino-cooled 
accretion. The neutrino-driven mass loss rate Mi, (interpolated 
between eqs. 1241 and 1261 ) is shown at the inner disk radius 
(r, = 10® cm; dotted line) and at the outer disk radius (near 
r^j, short- dashed line). The disk is advective to the right of the 
vertical line (eq. 1201 ). at which point the mass-loss will no longer 
be dominated by neutrino irradiation. 



as expected since Mi, oc r^^^. The vertical dot-dashed line 
marks where the disk transitions to being thick (eq. [20]), 
after which neutrino-heating no longer dominates the wind 
mass loss. 

Outflows that are launched from small radii, near r*, 
have the greatest potential to produce relativistic jets and 
to power high energy emission. But as we now argue, these 
neutrino-driven winds are too massive to become highly 
relativistic. Our calculation above focused on purely ther- 
mal, neutrino-driven winds, which accelerate matter to only 
a fraction of the escape speed (and thus are mildly rela- 
tivistic). However, in the presence of a strong, large scale 
open poloidal magnetic field, a more powerful, magnetically- 
driven outflow is possible. Magnetocentrifugal support in the 
wind's hydrostatic atmosphere may further enhance mass 
loss (e.g., Levinson 2006), but equation (|24p still represents 
the minimum mass loading on field lines which thread a 
neutrino-cooled disk. Figure [TOl shows that M^(rt) ~ 10~* — 
lO"'^ Mq s~^ during the thin disk phase. The luminosities of 
the prompt emission and late-time X-ray flares from short 
GRBs, however, do not typically exceed Lj ~ 10^° erg 

(and are often much lower; Nakar 2007). Thus, even 
assuming a modest radiative efficiency for the outfiow of 
Em ~ 0.1, the Lorentz factor F of a neutrino-heated disk 
wind must obey F ~ I/T,/[emA/^(r,)c^] < 5, which is incon- 
sistent with existing compactness constraints on short GRBs 
(Nakar 2007). A more likely source for the relativistic out- 
flows that power short GRBs and their late-time flares are 
nearly baryon-free field lines which thread the BH's event 



horizon (e.g., McKinney 2005). In addition, in ^5.21 we ar- 
gue that when the disk becomes advection dominated and 
neutrino irradiation effectively ceases, jet production may 
be more likely. 

5.2 Radiatively-Inefficient Thick Disk Winds 

At late times {t ~ fthick; eq. [20]) the disk transitions from 
thin and neutrino-cooled to being advective. At this point 
a neutrino-driven outfiow is unlikely to dominate the mass 
loss, in part because the neutrino luminosity precipitously 
drops (Fig.[3]&[4ll. In addition, because RIAFs possess a pos- 
itive Bernoulli parameter, a powerful viscously-driven out- 
flow is likely (Blandford & Begelman 1999; Stone & Pringle 
2001; Proga & Begelman 2003). 

In i|4.1l we showed that the disk becomes radiatively 
inefficient external to an "ignition radius" rign (x M^^^ 
(eq. [H]). The outer disk, near ra, thickens first (when 
rd ~ rign at f ~ tthick) and radiatively inefficient conditions 
move inwards as Md decreases. In the simplest picture, one 
might expect that the innermost radii become an RIAF only 
once Md drops from its value at t ~ tthick by an additional 
factor ~ (r^/rd)^^^ ■ In fact, the entire disk probably be- 
come radiatively inefficient on a timescale similar to fthick 
if the accretion rate which reaches small radii abruptly de- 
creases once the outer disk thickens (Fig. [BJ. Hence, at a 
time tthick, a significant portion of the accreting matter may 
be redirected into an outflow, with only a fraction ~ (r*/r£j) 
reaching small radii and accreting onto the BH (Stone & 
Pringle 2001). 

X-ray binaries typically produce radio jets upon transi- 
tioning from their "high-soft" (radiatively efficient) to "low- 
hard" (radiatively inefficient) states (e.g., Remillard & Mc- 
Clintock 2006). In analogy, once the inner disk becomes an 
RIAF, conditions seem to favor the production of relativistic 
jets (see also Lazzati et al. 2008) [f] 

Even if only a fraction (r*/rd) of the mass remaining 
when the disk thickens actually reaches the origin, the total 
energy supply available would be 

Ejct = ejctM(i(fthick)c^ I — 77- — r I 
\rd(tthick) J 

Q 1 nSO / Ejet \ / \ 6/17,, -4/17 / J'49\*'^"^^ 

~ 3 X 10 — — — T a' M,, ' ergs, 

VO.l/ VlOScmy 01 •J \^ 2 J ^ ' 

(29) 

where ejet is the fraction of the accretion energy used to 
power a jet and we have estimated M£j(tthick) and rd(tthick) 
using the self-similar thick disk solutions (eqs. |B5| and |B7] . 
respectively). Equation (|29|) shows that the accretion energy 
available from near r« following the RIAF transition is more 
than sufficient to power the late-time X-ray fiares observed 
following some short GRBs. If this is the case, fthick sets 
a characteristic timescale for late-time central engine activ- 
ity. If Q < 10~^, tthick niay be large enough to explain the 
~ 30 s delay until fiaring observed for some short GRBs (e.g., 
Berger et al. 2005; Villasenor et al. 2005). However, very late 
time energy injection, such as the Chandra flare observed 

* This is in stark contrast to jets powered by neutrino annihila- 
tion along the polar axis, which require a high radiative efficiency. 
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two weeks after GRB050709 (Fox at al. 2005), appears to 
require an alternative explanation. In addition, given obser- 
vational evidence for a ~ 0.1 in a number of environments 
(King et al. 2007), it may be more natural to associate iJjct 
and tthick with the energy and duration, respectively, of the 
short GRB itself, rather than the late-time central engine 
activity (see i|6]). 



5.3 Outflow Nuclear Composition 

The outflow nuclear composition has important conse- 
quences for the observable signature of compact object merg- 
ers. Nonrelativistic outflows are sufBciently dense to synthe- 
size heavy isotopes (Pruet et al. 2004; Surman et al. 2006), 
which may power transient emission via radioactive decay. 
The isotopic yield depends on the speed, thermodynamic 
properties, and the asymptotic electron fraction in the 
outflow[f| Although relativistic winds from the inner disk are 
unlikely to synthesize anything heavier than He (Lemoine 
2002; Beloborodov 2003a), is important in this case as 
well. A neutron-rich outflow may alter the jet's dynamics 
and the prompt and afterglow emission from that of the 
standard GRB fireball model (e.g., Derishev et al. 1999; Be- 
loborodov 2003b; Rossi et al. 2006). 

Figure [TT] delineates different regimes of outflow prop- 
erties and composition (as given by Y^) as a function of 
the wind launching radius r and accretion rate Md- We fix 
a — 0.1 and Al = 3Mq. The time-dependent evolution of 
the ring radius is shown for a solution with J49 = 2 and 
Md,o ~ 0.3Mq {solid line). At each time a given steady- 
state disk profile can be read off of this plot as a horizontal 
line that extends from the far left and ends on ra- Therefore, 
outflows from radii interior to rd contribute to the disk's to- 
tal nucleosynthetic yield. 

The ignition radius rign (eq. [19]) is shown in Figure [TT] 
with a short dashed line. For r > rign the disk is an RIAF 
and marked in the figure as "Thick Disk." In this case, a vis- 
cously driven outflow dominates ( i]5.2|l . Since outflows from 
RIAFs escape the disk in roughly the accretion timescale, 
these winds retain the midplane electron fraction (M08b), 
so that Ye" ~ Ve™"* <C 0.5 (because the disk itself freezes-out 
neutron-rich, as summarized in tj4.2l and Fig. [5J| . 

For r < rign, the disk is efficiently neutrino-cooled and 
marked in Figure[TT]as "Thin Disk." The absorption of neu- 
trinos, which heats the outflow and unbinds it from the BH 
may also alter its nucleonic composition. This drives Y^ to a 
value set by the neutrino radiation field Y^ , which in general 
is different from Ve™"^. A simple criterion was discussed by 
M08b for determining when Y^ ~ V/'. A typical nucleon 
in the accretion disk at radius r must absorb an energy 
~ GMruM /"ir to become unbound from the BH, so that 
Ni. ~ GMruM /2r{ev) neutrinos must be absorbed per nu- 
cleon. If we take A^i^ > Q ~ 2 — 3, then a typical nucleon 
has changed its identity (p ^ n or n ^ p) at least several 
times. 

This implies that all purely neutrino-driven outflows 
from radii smaller than 



^ The asymptotic electron fraction is germane because heavy nu- 
clei primarily form after freeze-out from /3-equilibrium. 
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Figure 11. Asymptotic electron fraction for disk winds as 
a function of the wind launching radius r and accretion rate 
Md (for a = 0.1 and M = 3Mq). The solid line indicates 
the location of the ring radius for our fiducial solution with 
Mdfi = O.3M0 and J49 = 2. The short dashed line is the "ig- 
nition" radius rign (eq. 1191 ). Exterior to this (marked "Thick 
Disk") the disk is advective with a viscously driven wind of com- 
position ~ y™''* < 0.5. Interior to rign (marked "Thin Disk") 
a neutrino-driven wind occurs. The dotted line shows r = 
with Q = 2 (cq. 1301 ) and determines where the neutrino ab- 
sorptions necessary to unbind matter alter the wind composition, 
so that ~ Y™'<i < 0.5 {Y^ ~ Y^) exterior (interior) to r^. 
The Md above which T[rj,) > 1 is plotted for BH spins of a = 
and a = 0.9. Above this line, the Pe and spectra differ and 
Y^ < 0.5, while below this their spectra are similar and Y^ ~ 0.5. 
In the region where r < Hgn, Ty(r*) < 1, and r < (i.e., the 
middle/lower left-hand trapezoid), y^f ~ Y^ ~ 0.5; these condi- 
tions are favorable for ^^Ni production (see SI5.4I I. 



GMm„ 



10'M3(eio>~'(Q/2)-'cm, 



(30) 



" 2Q(e,.) 

where (e^) = lO(eio) MeV, achieve Y^ ~ Y^ , independent 
of the disk's midplane composition. 

We plot with Q = 2 as a dotted line in Figure 1111 
where (ei,) is calculated from Md using our steady-state disk 
solutions (see i]5.ip . For r <^ r^, any neutrino-driven out- 
flow enters equilibrium with the neutrino radiation field (i.e., 
Yg" ~ Y"). For r > the outflow approximately retains the 
midplane electron fraction (i.e.. Ye" — Ye™"^). 

Although we have established the conditions under 
which Ye" is determined by neutrino absorptions, we must 
now address what sets Y^ itself. If the rate of neutrino ab- 
sorptions exceeds the rate of degenerate pair captures before 
the wind falls out of /3-equilibrium, Y^ is 



Ye"^ = 1 + 



L^, (e^,) -2A + 1.2AV(6^,) 
L,, (e,,) + 2A + 1.2A2/(6,,) 



(31) 



where A = 1.293 MeV is the neutron-proton mass difference, 
and Lv^/Lii^ and {^v^) / {e^^) are the mean v^jv^ luminosi- 
ties and energies, respectively, from a centrally-concentrated 
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source (Qian et al. 1993; QW96). Equation (|31|l demon- 
strates that the and spectra are crucial for setting 

^ e • 

Since the disk's luminosity and temperature peak at 
just a few Tg, is primarily determined by conditions at 
small radii. At early times, the accretion disk may be opti- 
cally thick near r* and so the Ve and i^e spectra depend on 
the temperatures at and i^e neutrinospheres, respectively. 
Since there are more neutrons than protons in the disk, the 
optical depth to fe through the disk is higher than to Ve', 
thus, the temperature at the z/g neutrinosphere is higher 
than at the i/e neutrinosphere. This implies La^ 3> L^^, 
(tiie) S> {^v^}, and thus Fe" ^ 0.5. Using 3-dimensional cal- 
culations of the merger of NSs with zero spin, Rosswog & 
Liebendorfer (2003) find that at ~ 15 ms following merger, 
L^, ~ 3.51/^,, (e^J ~ 9 MeV, and {e^J ~ 15 MeV, which 
implies Y" ~ 0.21, consistent with our arguments (see also 
Surman et al. 2008). We conclude that when the disk is 
optically-thick near r,, a neutron-rich outflow is again the 
most likely outcome. The critical accretion rate at which 
Tu{r,) = 1 is shown in Figure [TTI with a long dashed line for 
both a = and a — 0.9. 

Once the disk becomes optically thin near r, , the differ- 
ence between the i/^ and z7e spectra is much less pronounced. 
This occurs because (1) the neutrinos and antineutrinos 
originate from regions with the same temperature; (2) any 
net lepton flux out of the disk must remain modest (i.e., 
Lv^/{e^^) ~ Li)^/ {ea^)); and (3) the difference between the 
e~ and capture cross sections for kT 3> A — nieC? is 
small. Taking {e^^) ^ i^v^} 2> A, equation (|3ip shows that 
Y^ ~ 0.5, a value in the range required to produce ^^Ni 
(which we discuss further in tj5.4|l . Indeed, M08b used the 
steady-state, optically-thin Q-disk calculations of Chen & 
Beloborodov (2007; hereafter CB07) to calculate the neu- 
trino radiation fields carefully, and showed that Y" > 0.5 
over the majority of the disk (see their Fig. 1). Although 
the precise spectra extracted from an a-disk calculation 
should be taken with caution, the conclusion that the Ve 
and Ve spectra are similar for optically thin accretion (and 
Y" ~ 0.5) is probably robust. 

Figure[TT]illustrates that under most conditions the out- 
flows from hyper-accreting disks are neutron-rich. Neutron- 
rich material ejected during the initial dynamical phase of 
compact object mergers has long been considered a promis- 
ing source for producing Galactic r-process elements, whose 
precise astrophysical origin remains uncertain (Lattimer & 
Schramm 1974; see, however, Qian 2000). In addition, Sur- 
man et al. (2008) find that winds driven from the remnant 
accretion disk at early times (when it is optically thick; up- 
per left quadrant of Fig. Ilip are sufficiently neutron-rich to 
produce successful r-process. The outflows driven from the 
advective disk at late times, however, are unlikely to pro- 
duce r-process elements, given their modest entropies and 
electron fractions of Ye ^ 0.3 (Figs. [8] and |9]). Instead, this 
modest Y^ material will be synthesized to form intermediate 
mass neutron rich isotopes (Hartmann et al. 1985). 

5.4 ^^Ni Production and Optical Transients 

As summarized in Figure [TTI most of the material in the out- 
flow driven from a hyper-accreting disk will be neutron-rich. 
Nonrelativistic neutron-rich ejecta are difficult to detect be- 




Figure 12. Contours of total ^°Ni mass Mni = 
(XNi/0.4)_A/y^— 0,5 (in units of Mq) produced in the neutrino- 
driven outflows as a function of the initial disk mass q ^i^d 
initial ring radius r^fi, where My^=o.5 is the total mass loss with 
— 0.5 (based on the arguments in Fig. Illl l and Xni is the 
average ^^Ni mass fraction synthesized in the wind. The upper 
and lower panels correspond to non-rotating (a = 0) and rapidly 
spinning (a = 0.9) BHs, respectively. 

cause isotopes synthesized from low Ye material are them- 
selves very neutron-rich and typically possess very short 
half-lives, on the order of seconds (e.g., Freiburghaus et 
al. 1999). Thus, most of the radioactive energy is released 
at high optical depths and suffers severe adiabatic losses be- 
fore the photons can diffusively escape. By contrast, ejecta 
with Ye — 0.5 are easier to detect because they can pro- 
duce a significant quantity of ^^Ni (Hartmann et al. 1985), 
an isotope better suited to powering observable emission be- 
cause its half-life ~ 6 days is comparable to the timescale 
on which the outflow becomes optically thin. From Figure 
111! we see that outflows in a modest range of parameter 
space (middle/lower- left trapezoid) are capable of synthe- 
sizing ^''Ni. One caveat to this conclusion is that it only 
applies if the winds are primarily neutrino driven. If the 
outflow is instead magnetocentrifugally driven by a mod- 
erately strong open poloidal magnetic field (e.g., Levinson 
2006; Xie et al. 2007), then F^" < 0.5 can result, even if 
Ye — 0.5 (M08b). In what follows we assume that the wind's 
are primarily neutrino driven. 

Under this assumption. Figure [T2l shows the total ^^Ni 
mass, Mni ~ (XNi/0.4)My^=o.5, produced in outflows from 
hyper-accreting disks as a function of the disk's initial mass 
Md,Q and radius rd,o, where My^=o.5 is the total mass loss 
with Ye — 0.5 and Xni is the average ^^Ni mass fraction 
synthesized in the wind. We calculate My^^o.s by integrating 
the neutrino-driven mass loss (eqs. [24] and [26]) across the 
Ye — 0.5 region in Figure [TTI using rd{t) and Md{t) from 
the disk evolution calculations described in ij4] 

Pruet et al. (2004) present calculations of Xni which 
are parameterized in terms of the asymptotic entropy S" , 
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Figure 13. Luminosity of Ni docay-powered "macronovae" as 
a function of time since merger for Ni mass AfNi = 10-3 Mq 
and ejecta velocity v'^ =0.1 c. Light curves are shown for three 
values of the total ejected mass Mtot = 10" {solid line), 10~^ 
{dotted line), and IQ-^Mq {dashed line). The luminosities in V 
and J-Band (0.44 and 1.26 fira, respectively) are shown with thick 
and thin lines, respectively. The V-band upper limit on emission 
following GRB050509B from Hjorth ct al. (2005) is shown with a 
filled triangle. 

mass loss rate Aim, and asymptotic velocity n" of an out- 
flow with Fjf ~ 0.5f. is dominated by outflows 
from radii ~ 3 x 10® - lO'^ cm when Md ~ 0.1 - IMq 
s"^ (corresponding to I/52 ~ few); equation (|25|l thus gives 
~ 10 - 30fcs baryon"^ for the ejecta with Y^" ~ 0.5. 
Purely neutrino-driven winds achieve asymptotic velocities 
which are typically below the escape speed of the central 
object (e.g., Thompson et al. 2001); thus, the asymptotic 
kinetic energy is most likely dominated by energy released 
during the formation of heavy elements. Because ~ 8 MeV 
baryon"^ is released in producing Fe-peak elements, we es- 
timate that ~ 0.1 — 0.15 c. Applying these wind param- 
eters to Figure 3 of Pruet et al. (2004), we estimate that 
Xni ~ 0.2 — 0.5, thereby justifying our scaling for Xni in 
Figure [H 

Figure [T2] shows that for large initial disk masses 
{Md,o ^ O.IMq), the ejected Nimass, ~ 3x IO'^-IO'^Mq, 
can be appreciable. Disks with moderate initial radii rd,o ~ 
10'' cm are optimal for producing ^®Ni because they are 
sufficiently large to contain the radius ~ 10^ cm and 
yet are sufficiently compact to have a large initial accretion 
rate, which maximizes the neutrino luminosity and thus the 
neutrino-driven mass loss. Conveniently, initial disk param- 
eters from many compact object merger simulations (see 32]) 
are in the range required to produce ~ 10^* — 10~^Mq of 
Ni. 

The decay of Mm ~ 10"* - IO^^Mq can reheat 
the (adiabatically cooled) ejecta sufficiently to produce de- 
tectable transient emission. In order to explore this possi- 
bility, we calculate the light curves of ejecta heated by Ni 
decay ( "macronovae" ) using the method of Kulkarni (2005) . 
This simplified one-zone model accounts for the fraction of 
the gamma-rays produced by the Ni decay which are ab- 
sorbed by the expanding material (Colgate et al. 1980) and 



assumes blackbody emission at the photosphere, neglecting 
Comptonization. 

Figure [13] shows the V and J-band luminosities as a 
function of time since the merger for an outflow with Ni 
mass Mni = 10~^Mq which is expanding at v°' = 0.1 c. 
The V-band light curve peaks earlier because the tempera- 
ture at the photosphere decreases as the material expands. 
Somewhat after the peak in the light curves, recombination 
will decrease the opacity well below that considered here; 
thus our calculations are not quantitatively reliable at these 
times. The total mass Mtot ejected during the merger event, 
most of it neutron rich, is likely to be significantly larger 
than Mni; this provides additional opacity for the Ni-rich 
material. To explore the effect of this additional material 
on the detectability of the Ni decay, the light curves in Fig- 
ure [13] are shown for three values of Mtot: 10~^Mq {solid 
line), lO"'^ Mq {dotted line), and 10~"'^Mq {dashed line). As 
Figure [13] shows, larger Mtot: (1) delays the time to peak 
emission (tpeak is roughly cx M^^^); (2) increases the total 
fluence of the event by trapping a higher fraction of the 
gamma-ray emission; and (3) increases the peak wavelength 
of the emission, pushing it into the near-IR for large Mtot- 
We conclude that long wavelength (A ^ fim) observations 
at f ~ 1 day are the most promising for the detection of a 
Ni decay-powered macronova. 

Hjorth et al. (2005) place an upper limit of My > 27.5 
at t — 3.9 days on any emission associated with the short 
GRB 050509B (redshift z ~ 0.22); we mark this constraint in 
Figure [T3I with an arrow. For Mtot = O.IMq this constrains 
the ejected Ni mass to be Mni < 10~^ Mq (see also Kulkarni 
2005). As Figure [12] illustrates, compact object mergers are 
very unlikely to produce this much Ni, so the absence of a 
detection thus far is unsurprising. 



6 CONCLUSION AND DISCUSSION 

We have calculated the time-dependent evolution of accre- 
tion disks formed from compact mergers, and the properties 
of their outfiows. Since most of the disk mass resides at 
large radii, we approximate the disk as a ring at a given 
radius and calculate the dynamics and composition of the 
ring as a function of time. This ring model is calibrated 
to correctly reproduces the Green's function solution for a 
viscously spreading ring with viscosity v oc r^^'^ (appropri- 
ate for a thick disk; see Appendix [X)l . With this simplified 
model, we have studied the full parameter space of rem- 
nant accretion disks (different initial masses, compositions, 
etc.) and can follow the viscous evolution for arbitrarily long 
timescales. 

The energetics of the ring at a given time can be de- 
scribed by one of three models: (1) optically thick to neutri- 
nos and advective, (2) optically thin to neutrinos and geo- 
metrically thin, and (3) optically thin to neutrinos and ad- 
vective. A massive, compact disk (with a short initial viscous 
time tvisc.o; eq. [5]) will exhibit all three of these accretion 
phases, evolving from (1) to (3) as a function of time (Figs.[T]- 
[4|. Less massive disks, on the other hand, only pass through 
phases (2) and (3), or even just (3). Note that these phases 
refer to the energetics of the disk near the outer radius. At 
a given time, the disk may also undergo similar transitions 
as a function of radius; e.g., a disk that is advective at large 
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radii will be neutrino cooled and geometrically thin inside 
the ignition radius rign (eq. |19)). 

Neutrino-driven winds during the early-time optically 
thick and neutrino-cooled (thin disk) phases unbind so 
much mass that field lines connected to the disk cannot 
produce sufficiently relativistic material to power short- 
duration GRBs ( H5.ll and Fig. llOp . An alternative source 
for the relativistic material needed to produce short GRBs 
are nearly baryon-free magnetic field lines that thread the 
BH's event horizon (e.g., McKinney 2005). In addition, when 
the inner disk becomes advective {Md ^ Q.m al'^ Mqs'^ 
for a = 0), conditions appear particularly suitable for the 
formation of relativistic jets (by analogy to X-ray binaries, 
which produce jets when making a similar transition; e.g., 
Remillard & McChntock 2006; see Lazatti et al. 2008 for a 
similar argument in the context of long-duration GRBs) . 

Once the disk has transitioned to a late-time advective 
phase (phase 3 above), the properties of the disk become 
well-described by self-similar solutions. Ignoring for the mo- 
ment outflows from the disk, these solutions are rd oc t^^^ , 
Md oc t~^^^ , and Md oc t~^^^ . Power-law variations in the 
disk properties are a generic feature of a viscously evolving 
disk that conserves total angular momentum. These scalings 
are not, however, likely to be applicable in practice because 
outflows during the advective phase unbind most of the re- 
maining material ( tj4.1l fc |5.2|I . Energy produced by fusion 
to He and heavier elements also contributes to driving an 
outflow (Figs. |3]&|4}. Such outflows remove a signiflcant 
fraction of the angular momentum of the disk. This leads to 
a much more rapid decrease in the disk mass and accretion 
rate at late times (Appendix B3 and Fig. ^ . Signiflcant ac- 
cretion onto the central black hole will thus only last for a 
few viscous times after the onset of the advective phase. 

At the outer edge of the disk, the transition from a 
neutrino-cooled thin disk to the late-time advective phase 
occurs at a time tthick ~ 0.1 ^''/^^ (749/2)^^^^ s (eq. [20]). 
The rapid decrease in Md after the onset of the advective 
phase implies that the inner disk becomes advective at a sim- 
ilar time f H5.2l and Fig. |6}. Quantitatively, we flnd that for 
powerful winds with p = 1 (see eq. |B8) ). the inner disk be- 
comes advective at t ~ 0.2, 5, and 100 sec, for a = 0.1, 0.01, 
and 0.001, respectively (for our fiducial model with an initial 
mass of Q.IA/0 and an initial radius of ~ 3 x 10^ cm). Thus, 
for a ^ 10"'^, the timescale for the inner disk to become 
advective is comparable to the onset of observed fiaring at 
~ 30 sec in some short GRBs (e.g., Berger et al. 2005). Given 
the slow decline in disk mass with time before tthick, there is 
ample accretion energy available in the disk at this point to 
power the observed fiaring. However, there is observational 
evidence for q ~ 0.1 in a number of astrophysical disks (King 
et al. 2007); we thus doubt that tthick is large enough to co- 
incide with the onset of observed fiaring. Instead tthick is 
likely to be ~ 0.1 — 1 sec, comparable to the duration of the 
short GRB itself. In this case, the rapid decrease in the disk 
mass and accretion rate in the advective phase imply that 
the remnant accretion disk alone does not contain sufflcient 
meiss at ~ 30 sec to power the observed late-time activity 
from short GRBs, nor is there any physical reason to expect 
a sudden change in the disk or jet properties at this time. 

A more likely source of late-time flaring in compact ob- 
ject merger models is a continued inflow of mass at late 
times, such as is produced by the infalling tidal tail found 



in Lee & Ramirez-Ruiz's (2007) NS-NS merger simulations 
(see also Rosswog 2007). Similarly, the BH-NS merger sim- 
ulations of Faber et al. (2006a,b) show that ~ O.OSMq of 
material is ejected into highly eccentric orbits during the 
merger, which returns to the BH on a timescale > 1 s. How- 
ever, final conclusions regarding the quantity and ubiquity 
of late-time fall-back from NS-NS and BH-NS mergers must 
await fuU-GR simulations which include BH spin and real- 
istic EOSs. 

The second major focus of this paper has been on the 
composition of the accretion disk and its outfiows as a func- 
tion of time. For initial disk properties expected in com- 
pact object mergers (321), the disk typically comes into /3- 
equilibrium given the high temperatures and densities at 
small radii. As material spreads to larger radii, however, the 
composition of the disk freezes out before it becomes advec- 
tive at late times; at freeze-out the disk is modestly neutron 
rich, with an electron fraction Ye ~ 0.3 f i|4.2l and Fig. [9]). 
This neutron rich material - ~ 1O~^M0 for typical initial 
disk parameters - is blown away once the disk enters the 
advective phase at ~ tthick. These outfiows are particularly 
interesting given the low solar system abundance of mate- 
rial produced in nuclear statistical equilibrium at ~ 0.3 
(Hartmann et al. 1985). In a separate paper, we will study 
this nucleosynthesis and its implications in more detail. 

Although outfiows from compact object merger accre- 
tion disks are neutron rich in most circumstances, neutrino- 
driven winds from radii ~ 10^ — 10^ cm at accretion rates 
Md ~ 0.03 — 1Mqs~^ have electron fractions Ye ~ 0.5, 
precisely that required to synthesize significant amounts of 
^^Ni fFig. Ilip. We have calculated the total Ni mass ejected 
by compact object merger disks as a function of their ini- 
tial mass and radius ( >i5.4l and Fig. I12p . Disks with initial 
masses > O.IMq can produce up to ~ 10"^ Mq of ^^Ni. The 
radioactive decay of this Ni as the outflow expands to large 
radii will produce an optical and infrared transient peak- 
ing ~ 0.5 — 2 days after the merger, with a peak flux of 
i^L^ ~ 10^" ergs s"^ (Fig.[T3l). Because the Ni mass is likely 
to be a small fraction of the total mass of the ejecta (most 
of which is neutron rich), this transient is best detected at 
~ 1 ^m. As Figure [131 shows . current observational limits on 
SN-like transients coincident with short GRBs are about a 
factor of ~ 10 above our predictions. However, somewhat 
deeper limits from a moderately closer burst could start to 
put interesting constraints on short GRB progenitors. It is 
also possible that the decay of some neutron-rich isotopes 
could heat the outflow and contribute to the late-time ther- 
mal emission (although most such isotopes have very short 
half-lives). This possibility should be investigated in future 
calculations using a nuclear reaction network. 

Although we have focused on short GRBs throughout 
this paper, many of our results can be applied more broadly. 
For example, long duration GRBs show late-time activity 
and flaring similar to that seen in short GRBs (e.g., Falcone 
et al. 2007). For the reasons described above, this activ- 
ity is probably produced by a continued inflow of mass at 
late times (fallback from the stellar progenitor's envelope) 
rather than solely by the viscous evolution of the small- 
scale disk. As a final application of our results, we note that 
the accretion-induced collapse of a white dwarf to a neu- 
tron star (AIC) is expected to produce a compact disk of 
~ 0.1 — 0.5A'/q outside the newly formed neutron star's sur- 



© ???? RAS, MNRAS 000, [THIS] 



Time- dependent Models of Accretion Disks Formed from Compact Object Mergers 15 



face (Dessart et al. 2006). The calculations presented here 
describe the evolution of this remnant disk, with the one 
caveat that the composition of the disk in the AIC context 
may be strongly affected by neutrino irradiation from the 
newly-formed neutron star. 
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APPENDIX A: CALIBRATION OF THE RING 
MODEL 

The surface density E of an axisymmetric disk in a Keplerian 
potential with constant total angular momentum evolves ac- 
cording to a diffusion equation (e.g., Frank et al. 2002): 



dt 



3_9 

r dr 
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(Al) 



where v is the kinematic viscosity. Assuming that v depends 
only on radius as a power law, viz: u = Vfj{r j Rq)^ , equation 
(|A1|) is linear and, for an initial surface density distribution 
E(r, t = Q) — {Mo/2TrRo)S{r — Ro) which is narrowly peaked 
about the radius Ro, the solution (for n < 2) is given by 
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where Mo is the initial disk mass, x = r/Ro, r = t[12i'o(l — 
n/2)'^ / Rq], and 7m is a modified Bessel function of order m. 
For small argument y <^ 1, Im{y) takes the asymptotic form 
Im — (j//2)'"/r(m -f 1), where F is the Gamma function; 
thus, for late times or small radii such that r S> 2x^~"^'^ , 
equation (|A2|) reduces to 
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(A3) 



Most of the mass in the disk is located near the radius where 
the local mass Aid oc Er^ peaks; using equation (|A3|I , at late 
times this radius is found to be rpeak = -Ror^'''^""' . Hence, 
equation (|A3|) becomes valid near Tpcak for r 3> 1. 

The constant A, which relates the total disk mass at late 
times from the exact solution of equation (|A1|) to the mass 
defined by 7rE(rpeak)''pcak) can be calculated from equation 
~[3l) to be 



^(r > 1) 



Jo 



2n'Erdr 



7rE(rpcak)rp^ 



2e 



(A4) 



Similarly, the constant B, which relates the total disk angu- 
lar momentum at late times from the exact solution to that 
estimated by 7rErp^j,^(GMrpeak)^''^, is given by 



B{r » 1) 



7rE(rpcak)r 



5/2 
peak 



2e 



5~2n 



4-2n 



(A5) 



From mass continuity, the radial velocity is given by 
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which, using equation (IA3[l . gives the accretion rate at small 
radii 
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Equation (|A7|) is easily checked by noting that Alindt — 
Mo, which shows that the entire initial disk eventually ac- 
cretes onto the central object. In >i3.1l we introduced the 
following prescription for evolving the disk mass: 

fMd 



M, = 



(A8) 



where, in terms of the viscosity prescription adopted above, 
Uisc = r\/v = Uisc,o{rd/Ro)^~" and tvisc,o = Ro/^o is the 
initial viscous time. Assuming that the total disk angular 
momentum remains constant, J oc Mdr^^ = MoR^^ , the 



(A9) 



solution to equation (|A8|I is given by 

Md{t) = Mo[l + (4 - 2n)/(t/ivi=c,o)]"'/*''"'"' 



In our evolutionary calculations we set / so that the ac- 
cretion rate from the exact solution to equation (lAlf) {Min', 
eq. I A 71 1 matches the solution to equation (|A8|I at late times 
(i.e., in the self-similar limit). This requires 



/ = 3(1 - 7i/2)r[(5 - 2n)/(4 - 2n)Y 



(AlO) 



For an advection-dominated disk, v = aCsH oc Q,R^ oc 
r^/^; thus, n = 1/2, / ~ 1.602, A ~ 3.62, and B ~ 3.23. For 
a neutrino-cooled, optically-thin disk which is dominated 
by gas pressure, T oc r"^'^^" and v oc r^''^; thus, n = 6/5, 
/ ~ 1.01, A ~ 6.80, and B ~ 6.09. 

In Figure [Xn we show Mi-n/Md as a function of t/tvisc,o 
for n = 1/2 in order to compare the disk evolution derived 
from the exact solution of equation (lAlfl to that calculated 
from our simplified model. Figure lAll also shows the ratio 
of the total disk mass Mtot = /q 27rErdr calculated from 
equation HA2|) to the disk mass Md (eq. [A9j l of the sim- 
plified model, as well as the ratio of Tpeak (the radius where 
Er^ peaks, using eq. |A2j for E) to the radius determined by 
angular momentum conservation: rd = Ro{Md/Mo)^ ■ Figure 
lAll shows that, although the accretion rate in the two models 
differ at very early times (the initially narrowly-concentrated 
ring takes a short period of time to spread to small radii), 
they approach one another to <^ 20% by t ^ O.ltvisc.o. Like- 
wise, the disk mass and radii from the exact solution and 
simplified model are quite similar at all times. 

The numerical values for A and B given in equations 
(|A4|) and (|A5|) and employed in our calculations apply only 
to the mass and angular momentum distribution in the disk 
at late times (r 3> 1). Initially, the disk is entirely concen- 
trated at a single radius and A{t = 0) = B{t = 0) = 1; thus, 
A{t) and B{t) evolve significantly from early times until the 
disk enters the self-similar limit and so one might worry that 
the early-time description of the disk's evolution depends 
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Figure Al. Comparison of the accretion rate (solid), disk mass 
{short dashed), and disk radius (where the local disk mass 
peaks; dotted) as calculated from our simplified ring model to 
that derived from the exact solution of the diffusion equation 
for a (5— function initial mass distribution (eq. IA2n : we assume 
1/ oc r^/'^, as applies for a thick disk. The parameter / ~ 1.6 
(eq. lAlOh adopted in our model is chosen to ensure that the 
accretion rates match at late times (i.e., Mi^/M^ — » 1). Also 
shown is the ratio A{t)/B{t) (eqs. |A4) and |A5) ). a measure of 
the relative distribution of mass and angular momentum, which 
asymptotes to r[(5 — 2n)/(4 — 2n)] ~ 1.12 at late times. 



sensitively on the initial mass distribution. Our model only 
assumes, however, that the ratio A(t)/B(t) remains con- 
stant, which is a good approximation. To illustrate this, Fig- 
ure |Xl] shows A{t)/B{t) calculated from the exact solution 
(eq. [A2]) for n = 1/2. Note that A{t)/B{t) increases from 
unity to its asymptotic value A/B = T[{5 — 2n)/(4 — 2n)], 
which is ~ 1.12 for n = 1/2. 



Bl Neutrino-cooled, Thin-disk Solutions 

In the neutrino-cooled, thin-disk limit, the cooling is domi- 
nated by Urea, and the pressure is given by ideal gas. Com- 
bining local energy balance and continuity, Md — fAnvT,, 
allows us to solve for the temperature and column density 
as functions of radius. We substitute these into the angu- 



lar momentum equation, B{GMrd) 



1/2, 



jE — Jd, to solve 



for Md as a function of Md and Jd- We then assume the 
solutions have a self-similar form of Md oc f"*^, so that 
Md = -dMd/dt = jSMd/t. In this way we solve for = 5/8, 
Md{t), and subsequently any other variable of interest. The 
results are 



Md 



1.3 X lO'ViTs''* 



(Bl) 



Md 



2.7 X 10-7i".6 



5/8 



1-63.2 7 



-3/4 



m:. 



-1/4 



if) 



-13/i 



and 

Td = 4.1 X W'f!i'al^,^M-'^h'/\ 



Mqs'^ (B2) 



(B3) 



where /i.e = //1. 6, A3.6 = A/3.6, B3.2 = B/3.2, and t is 
measured in seconds, and the prefactors have been scaled 
to match our numerical results. The first thing to notice is 
that both Md and Md are rather insensitive to the choice of 
/ as long as it is near unity, and A and B only appear as 
a ratio, which is also nearly unity. This provides confidence 
in using this parameterization, and these specific values for 
the corresponding parameters, when the disk is not well- 
described by n = 1/2. This analysis also demonstrates the 
relative dependence on a. In Figure IBll we compare these 
scaling {dotted lines) with the numerical calculations. This 
shows that these solutions are only applicable for a short 
time. At times when t < tvisc the evolution is much flatter 
and is dominated by initial conditions. At later times the 
disk becomes advective and the solutions of the next section 
apply. 



APPENDIX B: 
SOLUTIONS 



ANALYTIC SELF-SIMILAR 



The late-time evolution of our disk calculations asymptote 
to power laws that are well approximated by analytic self- 
similar solutions. We derive these here to aid in interpret- 
ing our numerical results. Presentation is divided between 
neutrino-cooled, thin-disk solutions and late-time advective 
solutions. One could just as well derive analogous results 
for disks that are optically thick to neutrinos. We forgo this 
here since the initial viscous time is always sufficiently long 
that these solutions are never applicable to our numerical 
results. We conclude by presenting self-similar solutions for 
advective disks with substantial mass loss, since these differ 
significantly from the solutions without mass loss. 



B2 Late-time Advective Solutions 

In this limit, self-similar solutions can be found in an analo- 
gous way. The viscous energy release is carried by advection 
with the internal energy dominated by relativistic particles, 
so that 



Vr —aT . 

r 



(B4) 



SfAiv 

Combining this with mass continuity, gives the column 
depth as a function of radius, S(r) = {16/9ATva){M /r^Q). 
We then use this relation with B{GMrdY^'^Md = Jd and 
Md ~ PMd/t, to find /3 = 1/3 and the self-similar solutions 



Md 



3.7 X 10" 



\ 2 J 



A/q,(B5) 



Md = 1.2x10" 
and 



l(A3.(i\ -1/ 

(ff) -o:y^M3-/^ (if) t-/^Mos-\ (B6) 
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Figure Bl. Comparison of the numerical disk solutions [solid 
lines) with the analytic solutions for the thin, neutrino-cooled 
(dotted lines) and thick, advective limits {dashed lines). The nu- 
merical solution is the 0.3Mq disk from Fig. [T] 



ra = 2.3 x loVo^f Mg^^t^/^c 



(B7) 



These advective results are even more insensitive to A, B, 
and / than the thin-disk results. Equation (IB5|l - (|B7|l are 
plotted in Figure [Bll as dashed lines. The numerical calcula- 
tions follow these solutions very closely for times later than 
tthick (given by eq. [20]). 

Equations (|B5[) - (|B7|) can also be derived ignoring equa- 
tion (|B4[) . but assuming that the scaleheight is fixed at 
H/r ~ 0.6. This introduces the additional dependencies 
Md oc (i//r)-2/3, Md (X (H/r) -2/3, and ra oc (///r)'*/^ 
but gives nearly identical prefactors. 



B3 Advective Solutions with Mass Loss 



In H5.2\ we described how advective disks are likely to lose 
a substantial fraction of their mass to viscously driven out- 
flows. Because the outflow removes angular momentum as 
well - at least the specific angular momentum of the mass 
that is lost - the disk need not expand as rapidly to large 
radii. In addition, the disk mass and accretion rate decrease 
much more rapidly at late times than in the self-similar solu- 
tions described in the previous subsection. To quantify this 
effect, we follow Blandford & Begelman (1999) and assume 
that only a fraction ~ {r^/vd)^ of the available material is 
accreted onto the central BH. The remainder is lost to an 
outfiow. Thus the outfiow rate at any time is given by 

fMd 



= 1- 



rd 



(B8) 



We further assume that the angular momentum loss rate 
from the disk is given by 

{GMrdf^^. (B9) 



J 



-CMou 



where C is a constant that depends on the torque exerted 
by the outfiowing mass on the remaining disk. If the outfiow 
produces no net torque, an assumption that appears at least 
qualitatively consistent with the relatively small-scale mag- 
netic fields seen in global MHD disk simulations (e.g.. Stone 
& Pringle 2001), then the angular momentum loss is only 
that due to the specific angular momentum of the outfiow, 
and (Kumar, Narayan, & Johnson 2008) 

We solve equations (O, 0, (|B8[I . and (|B9|) . assuming 
A/B = 1 and y oc r^^^ (as appropriate for a thick disk). 
The solution depends on the relative magnitude of 1 — C 
and C(r^./rdY. For C{r,/rd)^ < 1 - C, which is true at 
nearly all times if equation ()B10|) is applicable, then 



rd ^ Tdfi 



l + 3/(l-C) 



t 



^visc.O 



2/3 



Md ~ Mdfi 



l + 3/(l-C) 



t 



^visc.O 



-1/[3(1-C)1 



(Bll) 



(B12) 



and 

Min 



r, 

rd,o 



l + 3/(l-C) 



t 



^visc.O 



-[l+3(l+2p/3){l-C)l/[3(l-C)] 



(B13) 



Note that if p = C = (i.e., no mass or angular momentum 
loss), then these self-similar solutions reduce to those of the 
previous subsection. However, for the case p = 1 consistent 
with a number of global advective disk simulations (e.g., 
Hawley & Balbus 2002), and in the absence of a net torque 
on the disk, C = 2/3 and these solutions correspond to 



Td oc t^l'^, Md oc and Min oc t~ 



8/3 



(see also Fig. [6]). 



This shows that the disk mass and accretion rate decrease 
subsantially more rapidly in time than in the absence of an 
outfiow, while the disk expands outward at roughly the same 
rate. If there is a net torque on the disk such that C ~ 1, 
then equations HBlHl - (|B13|) are not applicable. Instead, for 
C{r*/rd)^ S> 1 — C, the solution is given by (for p 7^ and 

t ^ iviac,o) 



rd(t)~ [(3 + 2p)M^^;g; 



1/(1. 5+p) 



t 



1/(1. 5+p) 



and 



Md{t) ~ Md,oexp[-Z?(t/tvisc,o)''''^'-'+^'] 



where 



D 



1.5 +p 



p(3-h2p)l-5/{1.5+p) ; I [/rPrl-g]l-5/(1.5+P) 



(B14) 



(B15) 



(B16) 



For p = 1 and for rdfi ~ r«, these solutions become rd(t) ~ 
r<i,o(t/tviBc,o)^'"' and'Md(t) ~ Md,o e^v[-^A^it/Uisc,o?^% 
The radius of the disk thus increases significantly more 
slowly, and the mass of the disk decreases much more 
rapidly, than in the self-similar solutions without mass-loss. 

The numerical solutions including mass-loss during the 
advective phase shown in §4.1 (Fig. [6]) assume that equa- 
tion (|B10|) is applicable and are indeed well-described by 
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the self-similar solutions given in equations (|BlHl - (|B13|l at 
late times. 
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